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Summary 

The  kinetic  theory  of  Prigogine  and  Balescu,  which  was  previously  developed 
for  a  homogeneous  plasma  in  a  magnetic  field,  is  now  applied  to  the  problem  of  dif¬ 
fusion  of  a  plasma  column  across  a  magnetic  field.  The  problem  of  deriving  a  kinetic 
equation  is  reduced  to  that  of  solving  a  two  dimensional  F 1  edholm-type  integral  equa¬ 
tion  with  a  complicated  kernel.  No  restriction  is  made  concerning  the  amplitude  or 
length  scale  of  the  inhomogeneity  across  the  magnetic  field.  The  results  are  preliminary 
in  nature,  and  suggestions  are  made  for  furthei  refinements.  Some  comparisons 
with  previously  existing  results  are  given. 
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SECTION  1 
INTRODUCTION 


The  phenomenon  of  diffusion  of  a  plasma  across  a  magnetic  field  is  of  great 
importance  in  the  problem  of  plasma  confinement.  There  have  been  many  theoretical 
papers  on  the  subject;  see  for  example  the  review  articles  by  Golant^  ^  and  Boeschoterf  l 
So  far,  however,  there  has  been  ve.y  little  work  done  on  the  kinetic  theory  of  inhomo¬ 
geneous  nonequilibrium  plasmas  in  magnetic  fields. 

We  report  here  some  pieliminary  results  on  the  diffusion  of  a  plasma  cr’  lmn 
across  a  uniform  magnetic  field.  It  is  believed  that  no  equivalent  results  have  yet 
been  published.  The  problem  of  deriving  an  equation  for  the  one  particle  distribution 
function,  including  a  proper  treatment  of  screening  effects,  has  been  reduced  to  that 
of  solving  an  integral  equation  of  the  F .  edholm  type,  with  a  complicated  kernel.  The 
method  used  is  an  extension  of  that  of  Prigogine,  Balescu,  and  coworkers.^^  ^  It  was 
applied  to  the  homogeneous  problem  in  preceding  works  ^  ^  ^  ^>to  be  referred  to 


as  papers  I,  II,  and  III  respectively,  We  rely  heavily  on  the  notation  summarized  in 

the  Appendix  of  paper  I  (  although  some  superscripts  J*.  and  z  are  suppressed,  and 

,  (4) 

x.  ®  is  replaced  by  Q).  Some  familiarity  with  paper  I  or  with  Balescu's  book  is 

J  J 

necessary  for  a  detailed  examination  of  this  report;  the  emphasis  is  on  the  differences 
between  the  present  calculation  and  that  of  paper  I. 

Although  the  plasma  is  assumed  to  be  homogeneous  along  the  magnetic  field,  an 
arbitrary  inhomogeneity  is  allowed  in  directions  perpendicular  to  the  field.  The 
infinite  order  perturbation  theory  is  tractable  for  this  case  because  the  guiding  centers 
are  constants  of  the  unperturbed  motion.  We  have,  however,  restricted  ou*  selves 
here  to  the  calculation  of  local  effects  in  a  single  species  plasma.  It  is  an  important 
problem  to  refine  this  model,  and  to  properly  take  into  account  internal  macroscopic 
electric  fields. 

More  details  concerning  the  physical  model  are  given  in  Section  2.  The  diagram 
perturbation  theory  is  summaiized  in  Sections  3,  4,  and  5.  The  main  results  are  given 
in  Sections  6  and  7.  They  are  compared  with  existing  kinetic  equations  in  Section  8. 


||)HiHI)j|||HH)|HHtilllll4tt)IHmwilllltl1ltl)IWI!Wll  ,:ii[W;'lii;ijlllH!!i:ilHmil)Wtilllll‘‘*lllilllHHlli:;Mllitlt!i^liHi||illlH)iliiiiiiMHii' . tiimii . .  . . . . . .  iiiiiiitiiiiiiiimniiimiiiiii 
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SECTION  2 
PHYSICAL  MODEL 


The  theory  presented  here  is  incomplete  in  the  sense  that  a  very  primitive  physical 
model  is  used. 

We  consider  a  single-species  plasma  in  a  region  of  volume  Q  containing  a  uniform 
magnetic  field  B.  The  length  of  the  region  in  the  direction  of  B  is  A,  and  the  cross 
section  is  £  =  Q/A.  We  will  assume  that  the  plasma  is  homogeneous  in  the  direction 
of  B.  In  addition,  the  limit 


N  -  oo,  A  -*  oo,  N/A  constant, 

£  arbitrary  (2.1) 

(cf.  eq.  (I.  2J.))  will  be  taken.  Thus  end  effects  (  for  example,  "short  circuiting")  are 
not  considered.  However,  an  arbitrary  inhomogeneity  in  directions  parpendic. nhir  to  B 
will  be  allowed.  There  will  be  no  restrictions  on  the  amplitude  or  length  scatc  of  the 
inhomogeneities,  although  these  scales  may  effect  the  plausibility  of  some  of  the  other 
assumptions. 

A  neutralizing  background  charge  of  unspecified  density  variation  across  B  is 
assumed  to  be  present.  This  is  an  ill-defined  artifice  introduced  solely  for  the 
purpose  of  restricting  ourselves  to  the  calculation  of  effects  due  to  relatively  local 
interactions  among  the  plasma  particles.  The  long  range  interactions  and  boundary 
effects  which  the  perturbation  cxparsion  of  paper  I  cannot  take  into  account  are  assumed 
to  effectively  cancel  each  other-  we  again  make  assumption  (1.3.26).  We  use  the 
Hamiltonian  (I.  2.  2)  with  W  set  equal  to  zero. 

In  the  homogeneous  case  discussed  in  paper  I.  the  generalization  to  a  more  realistic 
multi- component  plasma  would  involve  only  the  inclus’on  of  additional  indices  and 
summation  signs  in  appropriate  places.  In  the  present  case,  however,  such  a  generaliza 
tion  may  be  a  less  trivial  problem.  It  would  involve  a  careful  treatment  of  the  transverse 
electric  fields  arising  from  the  charge  separation  due  to  different  diffusion  rates  of  the 
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various  plasma  components.  This  should  lead  to  the  development  of  a  theory  of 
ambipolar  diffusion.  While  such  a  treatment  is  important,  it  is  beyond  the  scope 
of  the  work  reported  here. 

We  will  also  assume  that  the  one  particle  distribution  function  is  initially 
"gyrotropic" ;  that  is,  isotropic  in  directions  perpendicular  to  B.  This  may  restrict 
the  types  of  diffusion  that  the  theory  can  predict.  However,  the  theorem  that  the 
initial  condition  is  preserved  in  time  is  much  stronger  than  the  corresponding 
result  for  the  homogeneous  plasma.  In  that  case,  once  it  became  apparent  that 
the  higher  or  ler  initial  correlations  would  not  appear  in  the  final  results,  then  there 
was  no  way  in  which  a  preferred  direction  perpendicular  to  B  could  appear.  In  the 
present  case,  we  always  have  the  preferred  directions  corresponding  to  the  spatial 
density  gradients. 
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SECTION  3 

FORMAL  SOLUTION  OF  THE 
LIOUVILLE  EQUATION 

Our  starting  point  is  the  same  as  in  Section  I.  3,  except  that  we  do  not  introduce 
periodic  boundary  conditions  with  respect  to  the  \  ®  X  }  (hereaftjr  called  ). 

The  fundamental  reason  why  the  treatment  of  arbitrary  inhomogeneities  perpendicular 
to  B  is  tractable  has  already  been  noted  in  the  paragraph  following  eq.  (I.  3. 13).  The 
^  Qj  J  are  constants  of  the  unperturbed  motion,  and  hence  can  be  treated  in  the 
same  way  as  the  velocity  cylindrical  components  ^a^  ,  p^  Z  ^  .  (Fro—i  now  on,  the 
superscript  z  will  be  suppressed.  )  We  had  no  reason  to  assume  any  sort  of  homogeneity 
in  velocity  space  (except  through  the  condition  of  staoility);  hence  none  is  needed  with 
respect  to 

Thus  eqs.  (I.  3. 14)  and  (I.  3,  18)  are  replaced  by 


yt)Hy)2j|,|Q|.iaj|,!ej|,|Pj|>t) 

=  I  Z  o'  (t)exp  [i£  (k.z.  -  n.9,)], 

|k|{„|  WM  J 


(3.1) 


1  oo 

‘  zm  i  dc  l  I  (-«  >’ 


q=o 


[(I  kj<0)pj/m)  +  “(Inj(°,)  '  C1 


(°h  rn-1 


I  l  l  <bj,e’,Hnj(0)l  l6Lhul  l^'11  Hnj(1)| 

|k(11|n<1'|  h<d 
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[<£  k.(1)p./m)  +»1(^n,|1))  -  C]'1' 


I  l  l 


ik<q)[  |»(q)|  i< 


m 


<  |k.(‘>-1)|  |n.lq-1)t  |6L.J  |k.(q)||n.(q)|> 


['Ikj  Vm,  +  “’lIni  ,'cr  Pik'q>H„«»l(0) 


(3.2) 


J 


J 


Here,  the  k's  are  restricted  to  values  of  the  form  (2n/A)(integer),  and  the  n's  are 
integers. 

The  matrix  elements  are  defined  as  follows: 


u.ii^i  i6Lhuiikj ! i ni i  > 


=  A-N(2n)-N  J dNz  dN0 

exp  [-i  l  (k.z.  -  „.9.»  6Lhu  exp[i  £  (k^  -  ^  6.)]. 

j  j 

They  are  operators  with  respect  to  ,  (a.\  ,  and  .  An  explicit 

J  J  J 

expression  analogous  to  eq.  (I.  3.  25)  can  be  calculated  in  the  same  way: 


=  C  n  SKr(k  -  k')  SKr(n  -  n')] 
j(*h,u)  33  33 
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OO  £TT 

6Kr(k,  +k  -  k/  -  k')  \Ja/\  dB 
'  h  u  h  u  Jo  .to 


[-rrmiuA  i/1  +  [kh  -  k' ]2  +  h2)]"1 


exp  [i(nh  +  nu  '  \  *  »'  )  8]  {  Jn  -  n'  (^h) 

*  h  h 


/.(^a  )  [muu(k  -  k')  (  -  T“) 

u)  u  h  h  *ph  aPu 


-(n  -  n  )  u  h  h  op,  op 


n,  n 

/  /,  _1_  *z  4.  Jl  _2_  +  Ji.  i 

(  -  X  9Q,  J  +  a  9a  +  a  *a  ^ 

=hu  h  h  u  u 


-J  ,  '1  (^a  )  r1  dr  J n  n 

-(n  -  n  )  u  a  3a  n  -  n,  « 

a  a  h  n  n  n 


\  '  \  (^h)  a- Sa-  J'(n- '  n-]  (^u)} 


u  “  u  u  u 


exp  [i  i  *  Qhu]. 


(3.3) 


where  Q  =  Q  -  Q  ,  and  ^ is  the  vector  defined  by  the  polar  coordinates  8 
—  hu  —h  — u  - 

(y^B  =  o). 

It  is  possible  to  perform  the  integration  over  B.  Although  the  result  will  not  be 
used,  we  write  it  for  completeness: 


<{^|  |n.  |  |6LhJ)k:|  |n'|  > 


=  [  n  6Kl  (k.  -  k.')  b*1  (n.  -  n')] 
j(4h,u)  J  J  J  J 
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1/  -CO 

6 1  (kh + ku  kh-k:)2nI 


[-nmu)A{«^2  +  [kh  -  k^]2  +  H2)]  1 


{ 


[mu,(k  -  k')(- —  -  - — )  -  —  - — 
h  h  ap,  *p  a,  9a 
h  u  h  h 


/  ,ii 

n  „  n,  +  n  -  n.  -  n  . 
u  h  u  h  u  o 

a  *a  Q,  *0, 

u  u  hu  hu 

J'  /  ,  Q,  ) 

n,  +  n  -  n,  -  n  hu 

i  /  h  u  h  u _ o 

Qu  J  ,  /  /  (/Q,  )  J 

hu  n,  +  n  -  n,  -  n  *'  Tiu  hu 
h  u  h  u 


[J  /  (jf  a,  )  J  ,  / «(  ^ a  )]  +  [J  /  ) 

nh  ‘  "h  h  ‘(nu  ’  nu)  U  "h  '  "h  H 


u  u 


h  h  u  u 


J  /  ,  (V^Q,  )  exp  f  i(n,  a  n  -  n'  -  n')($  i 

n  +  n  -  n,  -  n  hu  h  u  h  u  Thu 

h  u  h  u 


where  =  arctan  (Q^  /  Q^*  )  (cf.  eq.  (I.  3.8)  ff). 

(9) 

We  also  note  that,  according  to  Barabancnkov  ,  eqs.  (6)  and 
over  £  is  simple  if  >  a^  +  a  .  One  should  however  check  the 
those  formulae  before  applying  them. 


y)],  <3-4) 


(7),  the  integration 
range  of  validity  of 
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SECTION  4 
INITIAL  CONDITIONS 

For  the  physical  model  described  in  Section  2,  the  arguments  in  Sections 
I.  4  and  I.  5  go  through  almost  without  change.  All  that  is  necessary  is  to  replace 
fl  or  A/8tt3  by  A  or  A/2TT  in  various  places,  and  to  include  integrations  over  j  Qj  } 
with  the  velocity  integrals. 

The  reduced  distribution  functions  can  be  defined  as  before: 

f  r ,  s(zT  *  *  *  '  Zr’  I9a|  ’  iaa  I'  i9a  M  Po  M 


1! _  P; 

-  r)!  J„ 


N! 

(N  -  r) 


•N-  r  f'  dN’SQ  dN-SadN-S0dN-8p 


(1,.  .  •  ,  r) 


J. 


M 


N 


2  2 


i  n  .1  QjM  V'l  v  -t  pjl  '  t>; 

h=  s  +1  J 


(4.1) 


cp  =  f  ,cp  =  cp  =  f  ;  0  <  r  <  s  <  N; 
8  O,  S’  1  0,1  - 


with  the  normalization  condition 


J  dS  Q  dSa  dS9  dSp  (H  m  uu  a^)  cpg({  a  }  ,  t)  =  1. 
j  a 

They  can  be  expressed  in  terms  of  the  p  * s : 

frs  =  AN(N»'A)r  l  l  Vs 

k.  ..k  in  | 

1  r  *  at 


(4.2) 


p ,  I.. .  r  j  a  |  exp  (.j^  ^  0p) 


exp(i[kjZj  +  .  .  .  +  k^z^D, 


(4.3) 
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which  is  analogous  to  eq.  (I.  4.  7).  Here, 


f'  dN'5  Q  dN'Sa  dN'Sp 

>1 


(4.  4) 


The  Ursell  correlation  functions  Ug(|o  |  )  are  defined  by  eqs.  (1.4.11),  with 
Q  replaced  by  A. 

Our  object  is  to  obtain  equations  for  the  one  particle  distribution  function  cp. 
Initial  conditions  analogous  to  those  of  Section  I.  5  are  imposed: 


f 

l 


(0)  =  f 
s  r.r 


(o)  n 


+ 1 


^(QT»aT>  6t>px;°). 


(4.5) 


cp(Q, 


pt-°)  = 


:(9t 


(4.6) 


Ug(  {  o  }  ,  0)  is 


finite  in  the  limit  (2. 1), 


and 


lim  U  (.  .  .  i  0. ,  .  . 
s 

\  -*  00 


A 


.  .  ;0)  =  0  for  Ua 


- 1 


cor  r 


(4.7) 


It  is  assumed  that  there  are  no  initial  correlations  among  the  guiding  centers  or 
velocities  which  are  independent  of  the  distances  |  |  z.  -  |  |  ,  in  the  limit  A  oo. 

The  condition  (4.6)  is  stronger  than  for  the  homogeneous  case,  because 
of  the  preferred  directions  corresponding  to  density  gradients  across  B.  However, 
a  strong  theorem  can  be  established;  namely,  that  the  condition  is  preserved  in  time 

in  the  modified  weak  coupling  approximation. 

We  can  then  easily  establish  the  analogue  of  the  result  (I.  5.  22):  If 
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1.  .  .  r  jo}_  .N..  ,,„.r-w  ,l...r  jo} 

k**-q  {n\  k.  ..q  |n[ 


k+. .  ,  +q  =  0,  k  $  0,  . .  .  ,  q  ^  0, 


(4.  9) 


then 


(2tt/ A)W  V..  pf  '  *  *  r  |°|  (0)  remains  finite  in  the  limit  (2. 1),  (4.10) 

^JN-s  k.  ,.q  j  n  1 

where  w  and  w#  have  the  same  meaning  as  before.  This  follows  from  expressions 
such  as  that  analogous  to  (I.  5. 13): 


U2<1. 2;  0)  =  (21T/A)  2  l  Vz 
k(*0)  |  n  | 


pk  -k  n2  <°>  "a  e„  ]  [ik<zi  -  s2>  ] 


(4.11) 
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SECTION  5 

SELECTION  OF  THE  DOMINANT  TERMS 

The  diagram  representation  of  the  formal  solution  (3. 1),  (3.  2),  (3.  3)  is 

almost  identical  to  that  for  the  homogeneous  case.  The  "semi-connection"  property 

(I.  7,7)  remains  valid,  because  although  there  is  an  extra  term  in  (3.  3)  containing 

d/5Q  .  there  is  no  dependence  on  Q  or  Q  to  the  left  of  this  factor. 

“  —hu  ft  u 

By  calculating  the  A  dependences  of  the  various  diagram  contributions  to 

p  .  *  !  °  !  (t)  (using  condition  (4.10)),  we  can  establish  that  condition  (4.10) 

N-s  j  k  }  |  n  I 
remains  valid  for  t  >  0. 

J) 

As  for  the  e  dependence  of  the  initial  distributions,  there  is  no  essential 

r  n  g  n  r  n  n 

change.  In  equations  such  as  (I.  7.  20),  one  replaces  J  d  x*/0  by  J  d  z/A  ,  and 
(4tt^/E)V  by  J  (It  is  then  possible  to  integrate  over  8  to  obtain  a  simpler 

result;  cf.  eq.  (3.4))  We  therefoie  assume  that 


,  ...  ,  2. r-w-w  ,  ... 

|k}  {  n}(0)  “  (G  ]  P  1 0  ||  0|(0)’ 


(5.1) 


which  is  similar  to  eq.  (1.7.24). 

Let  us  wiite  the  equation  for  the  contribution  to  Cp(a,  t)  corresponding  to 
(I.  81)  or  (I.  8.  3): 


cp(a,  t)=...  -  (Z"*)’1  j>  dC  exp(-iCt) 


(-c''‘  III  <-Ne2> 

k  JDJD 

n0  "h 


t  0  H  0  I  l6Lhalkh=-k’ka=k’41)’41,-)°l'  > 


(k[pa  -  ph]/m  +  w[n(1)  +  n^]  -  C)_1  £  (-Ne2) 

"(2) 


u 
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<V-k'Vk’  "i1*'  4"  l°l'  |6Luh|ku 


*k-ka  =  k> 


■£’.  n«2,.t0|'>(K[pa-Pu3/m+»[n-+„-l 


ID  ^  C)-i 


(■e2,  <  ku  =  k,  ka  =  k,  „J»,  «»> ,  t  0  |  '  |6LaJ  1  0  H  0  I  > 


.'I 


('°  *  °|  0  I  |  0  }(0)’ 

where  {  0  |  '  means  that  all  suppressed  variables  are  zero.  We  replace  n^ 

(1)  n(21  k  n  .n  -n  ,  and  make  the  following  abbreviations  in  notation: 

h  ’  u  7  a  h  u 


(5.2) 


(1) 


<0  |ah  |  >  s  -  Ne2i  ^  £ 


n  n, 
a  h 


:  1 0  II  0  I  l6Hjkh=-k’  Vkl  na)  =  V  nh’=-V*  0t,> 

(with  certain  terms  that  vanish  under  omitted) 


oo 


2tt 


“I  I  I  A  d  I  d,h  exp  (i["n  '  na]  Sh> 


n  n,  o 
a  h 


2Nexp(i/  h-  Q.„a)  j 

- - - -  Lmyjk  - — 


iTTmuuA^^  +  k2  +  h2) 


SPn 


K  — — —  )z  +  —  —  ]J  (^,a  )J  (^  a  ), 
'  ~  h*  3Q  '  a  5a  J  n'  ha'  n'  hh 
—  a  a  a  u  n 


(5.3) 
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•[muuk( 


dp 

r  u  a 


)  -  ( 


x — —  )z 
—  ua 


_ u  _d_ 

a  da 
u  u 


_a _ d_ 

a  da^ 

a  a 


Then  eq.  (5.  2)  becomes: 


co(a,t)  =  ...  +  j>  dC  exp(-if.t) 

c 


OCV1  f-  )_j<0  M  ><i[Gah  ■  C])_1 


«■  M>  (iCG^u  -  CD'1  <  |ua|  0  >  Pj  o  ||  o  j  (0), 

where  P|  Q  j  |  Q j.O)  =  A N  P  ^  0  J  j  0}  (°)  can  be  replaced  by  n  cp(j,0)  (cf.  eq.  (1. 10.  3)). 

The  arguments  from  cq.  (1.8.  8)  to  eq.  (I.  8. 18)  can  then  be  repeated.  In  the 
asymptotic  limit  t  >  >  t  ,  eq.  (5.  6)  becomes 


(5.5) 


(5.6) 


CD 


co(a,t)  =  ...  +  t  J  dk  <  0|ah|>6.(Gah)<|hu|>6jGau)<|ua|0>ncp(j,O) 

-oo  j 


(5.7) 


In  selecting  the  dominant  terms,  we  again  make  the  modified  weak  coupling 

approximation.  Out  of  all  the  diagrams  proportional  to  particular  powers  of  N/A 

2 

and  t,  we  select  the  ones  proportional  to  the  lowest  power  of  e  .  Then  the  arguments 
of  Section  I.  9  can  be  repeated  (see,  however,  the  footnote  at  the  end  of  Section  H.  2). 
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The  dominant  diagrams  are  semiconnected  chains  of  ring  diagrams 
In  particular,  we  have  the  result  analogous  to  eq.  (I.  9.  9): 


If  k':x}  =  0, 


then  n. 
J 


(M 


0,  for  any  j  and  X. 


This  implies  chat  the  isotropy  condition  {4.  6)  is  preserved  in  time. 
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SECTION  6 

THE  SINGULAR  INTEGRAL  EQUATION 

The  major  mathematical  complication  arising  from  the  presence  of  the 
inhomogeneity  across  B  is  already  apparent  from  eq.  (5.  7).  The  matrix  elements 
contain  two  labels  instead  of  one,  corresponding  to  the  labels  to  the  right  and  left 
of  the  diagram  vertex.  This  feature  cannot  be  rectified  by  shifting  the  position  of 
the  Bessel  functions,  as  was  done  for  the  homogeneous  case,  because  of  the  presence 

of  different  dummy  variables  h*  ^u'  a°d  ^  instead  °f  a  single  variable  -* 

This  complication  seems  to  prevent  the  derivation  of  a  kinetic  equation 
of  the  same  simple  form  as  for  the  homogeneous  case.  However,  some  progress 
has  been  made,  as  will  be  seen  below. 

It  can  be  shown  that  the  analogue  of  eq.  (I.  10.  1)  is 

cp(a,  t)  =  .  .  .  +  ^n-1  1  i  dk  *  °^a1^  > 

oo 


•<  |  34 1  >  6JG42)<  1 45 1  >  6JG52)<  1 26 1  > 
•6_(G66)  <  1,57 1  >  6JG76)  <  168|>  6_(G?8) 

■  <  |79  |  >  6_(G98)  <  |  8  10  i  >  6JG9  1Q)  <  j  10  11 1  > 
•6JGq  u>  <  '.9  !2|  >  MG12  u)  <  |  11  13  |  > 

•6  (G.,  .,)<  1  13  12 1  0  >0  cp(j,0), 
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where  the  asterisks  denote  the  complex  conjugate.  The  arguments  leading  to  eq. 

;I. '  .15)  are  not  affected  by  the  more  complicated  expressions  for  the  vertices;  the 
cha  ->£  ring  diagrams  can  be  factorized  as  before.  In  particular  the  condition  (4.5) 
is  preserved  in  time,  at  least  for  r  =  0. 

However,  we  cannot  express  the  results  in  a  form  exactly  analogous  to  eqs. 
(1.10.18),  (1.10.25).  Instead,  we  define  a  new  quantity  ^(a,  y),  which  is  related  to 
the  binary  correlation  function  by  a  Fourier  transformati"  .1: 

oo  — i 

U_(a,y)=  j  dk  )  S  F(a,y)  exp(ik[z  -z  ])  exp  (-i  [n  fl  -  n  0  ]). 

2  J-oo  L  “  aY  a  a  Y  Y 

n  n 

a  y 


(6.2) 


(We  suppress  the  parameters  k  and  t  from  the  arguments  from  now  on,  since  the 
equations  no  longer  involve  the  initial  time  t  =  0.  )  Then  it  can  be  shown  that 


=J  dk  J  <0  |oth|  >  F(a,  h), 


-oo  h 


(6.3) 


F (a,  h)  =  q(a,  h)  +  6_(Gah)  J  <  ihul  >  ^(h)  F(a,  u) 


ft 

+6.(Gah)j  <laul  >  cp(ct)  F  (u,  h), 


(6.4) 


where 


q(a,  h)  H  6_{G  )  <  I  ha  I  0  >  cp(h)  cp(a), 


(6.5) 


OO  00  ..  _ 

r  r*  f>  p  2  2 

i  f  s  :  dQ.  I  da.  dp.  2tt  m  uu  a.f. 
Jj  JE  “J  o  ]J-oo  1  1 


(6.6) 
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Note  that  if  the  matrix  elements  did  not  depend  on  the  second  labels,  then  cne  could 

express  Scp(a)/dt  in  terms  of  [*  F(a,  h),  and  derive  an  integral  equation  for  f  F(a.  h). 

Jh  h 

The  singular  integral  equation  (6.4)  is  very  complicated,  and  docs  not  seem  readily 
amenable  to  the  treatment  described  in  Section  II.  4.  By  making  transformations 
similar  to  those  in  Section  7  below,  we  could  eliminate  some  of  the  variables  by  a 
sort  of  "barring"  operation,  and  derive  a  very  lengthy  equation  for  a  quantity 
$(a,  q^)  which  is  independent  of  a^.p^,  and  which  determines  F(a,  y)  through  the 
following  relations: 

F(a,  Y)  =  J  dqa  j  dq  expti^  ■  ^  •  Qy)]  F(a,  y)  , 


(6.7) 


and  simila  ly  for  qfe.,  y)  and  q(a,y); 


f,  .  A.  ,  ,  4  ne^N  .  ,r  .  r”  qedqc 

F(a’ Yl  ■  q(a’v» +  wr  6-'V  j  r-j~  i 

o  q  +  k  +  h 


J  Jn  KV  •=*p(-i,y£)[mj-k¥r 


+  i[aE*aY]z+r  -£r  ]  $(V!  !sy  ■  ae  -  iY» 
y  y 


)  J  ,  (q  a  )  J  (q  a  )  exp  (in  *  )  $(a,o  ) 
L  n  L  v  a  a  v  yx  a  a  a  -*e 

v  a  a  a 


+  a  < — >  y 


(6.8) 
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where  ci>  is  the  Fourier  transform  of  cp: 

cp(Y)  =  Jdky  exp  (iky  •  Qy)  $(Y),  (6>9) 

and  £e  is  defined  by  the  polar  coordinates  q  it  Such  transformations  have  been 

used  extensively  by  Rostoker  and  coworkers,  and  others;  see  for  example  Reference 

(10). 

We  do  not  bother  to  write  out  the  equation  for  *(d,a  ),  because  a  simpler 
set  of  equations  can  be  found  in  another  way.  In  the  next  section,  we  use  the  sum¬ 
mation  procedure  developed  by  Resibois  (U)  in  order  to  derive  equations  equivalent 
to  (6.  3),  (6.  4),  but  having  a  more  explicit  form. 
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SECTION  7 

THE  FREDHOLM- TYPE  INTEGRAL  EQUATIONS 

The  Resibois  summation  procedure^)  has  been  described  very  clearly  and 
concisely  by  Balescu  (Reference  (4),  Appendix  10).  Singular  integral  equations  are 
avoided,  and  kinetic  equations  are  obtained  diiectly.  The  method  is  largely  independent 
of  the  nature  of  the  matrix  elements  (vertices).  It  generalizes  immediately  to  the 
present  physical  model,  provided  that  we  can  consider  the  integrations  term  by  term 
in  the  n's.  The  n's  are  held  fixed  until  we  actually  sum  the  "primitive  diagrams". 

The  complication  mentioned  in  the  last  section  also  appears  here,  but  we  are 
able  to  proceed  much  farther.  We  do  not  get  simple  geometric  series,  and  therefore 
no  simple  kinetic  equation  is  obtained.  However,  we  can  reduce  the  problem  to  that 
of  solving  a  two  dimensional  integral  equation  of  the  Fredholm  type  but  with  a  complicated 
kernel.  Further  progress  has  not  been  seriously  attempted  at  the  time  of  writing;  it 
would  seem  that  some  approximation  may  be  necessary  for  the  derivation  of  a  kinetic 
equation  in  a  closed  form. 

We  replace  SJG^)  bY  the  equivalent  form  CifG^  -  iG)]’1,  and  write,  according 
to  the  Resibois  theorem, 


V*  h  f  .d£'  I 


I.s 


*<  |  SI |  0  >  TT  cp(j), 
J 


(7.1) 


where 


\i(C/)  =  [i(Ga  -  C'>]‘1<M>' 


*[i(G?  -  C')]'1  <  1 23  |  >*  ...  Ci(Gj  -  G')]  , 
1/hS(i€  -  C ' )  =  ti(-Gh-  ic  +  C')]'1  <|H  1+  1  |> 
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[i(-Gj  +  j-  ie  +  C')!  1  <  j I  +  1  I  +  2  |  > 

...  [i(-G  -  ic  +  C')]'1  • 
s 

Here,  C'  is  a  path  of  integration  which  consists  of  a  straight  lint  from  £'=-®+ie/2  to 
Q'  s  f  oo  +  ic/2,  together  with  the  usual  semicircle  in  the  lower  half  plane.  (The 
orientation  is  opposite  to  that  of  path  C  in  eq.  (5.  2).  ) 

From  eq.  (7.1),  it  is  a  straightforward  matter  to  obtain  recurrence  relations, 
and  the  resulting  integral  equations: 

F(a,h)  -  ~  q  dC'  7(a,h)  (7-2) 

c 

h)  =  [i(G&  -  '')]  1  [A (° .  h)  +  j*  <  juu  j  >  '■p(a)r(u,h)  ]  , 

u 


A(a,h)  [i(-Gji~i-:+  C,)1  1  rt(a,h)+  f  <  |hu  ]>cp(h)  A(a.  u)  ]  , 

u 

(7.4) 

8  (a,  h)  =  <  |hajo>  "f(h)rp(a).  (7.5) 

It  is  claimed  that  the  sot  of  equations  (6.  3),  (7.  2),  (7.  3),  (7.  4),  (7.  5)  is  simpler  than 
the  set  (6.  3),  (6.  4),  (6.  5).  Note  that  if  one  can  solve  eq.  (7.  4),  then  there  is  no 
trouble  in  solving  eq.  (7.  i),  because  it  has  essentially  the  same  kernel. 

Let  us  now  investigate  briefly  the  .  integration  in  cq.  (7.  2).  It  is  relat  vely 
easy  to  perform  the  integration  for  the  [i(Ga  -  C')]’  A  (a,  h)  contribution  (eq.  (7.  3)). 
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The  factor  [i(-Gh  -  ie  +  C')]"1  in  eq.  (7.  4)  has  no  poles  in  the  area  enclosed  by  Cf 
Hence  it  is  very  plausible  to  assume  that  A  (a,  h)  has  no  singularities  there  either,  at 
worst  there  is  some  sort  of  stability  condition  involved.  We  can  write 

s 

—  £  dC#Li(G„  -  C,')]'1  A(a,  h)  =  A(a,h)L,  _r 

2tt  J^/  a  ^  '  a 


s  A  (a,  h), 


where 

A  (a,  h)  =  6_(Gah)[0(a,h)  +  <|hu|>  cp(h)  A(a,u)]. 


(7.6) 

(7.7) 


The  other  contribution  to  L’(a,h)  is  more  difficult  to  calculate.  We  merely  note  that 
the  contribution  from  the  semicircular  portion  of  c'may  be  neglected,  and  then  change 
notation.  Replacing  Q1  by  4  ie  /  2 ,  we  get'. 

,  oo  .  ^ 

F  (a,  h)  =  A  (a,  h)  +  j  ^  dGx  C  6JGax)  Ju  <  I aU  I  >  r  ^u’  h:  ’ 

(7.8) 

r(a,h;x)  =  6_(G,  .)  [A(a,h;x)  +  f  <  |  au| >'°CD(a)  Hu,  h;  x)]  , 

(7.9) 

A(a,h;x)  =  6JGxh)  [  0  (a.h)  +  f  <  |hu|  >  o(h)  A(a.  H;x)  ]  , 


(7.10) 
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an 

A  (a,  h)  =  A  (a,  h;a)  =  J  dGx  6(Cax)A(a,h;x).  (7.11) 

-oo 

We  now  examine  the  problem  of  simplifying  the  integral  equations(7.  3)  and 

(7.4),  or  (7.9)  and  (7.10).  All  these  equations  have  essentially  the  same  kernel,  and 

can  in  principle  be  solved  in  the  same  way.  We  consider  eq.  (7. 10)  as  an  example, 

and  abbreviate  A  (a,  h;x)  by  A(h),  since  a  and  Gx  appear  only  as  parameters.  Also, 

6  (G  )  ©  (a,h)  will  be  written  as  -•  (h),  so  that 

—  xh 

A(h)  =  -(h)  +  6_(Gxh)  Ju<|hu|  >cc(h)  A(u).  (7.12) 

The  following  expansions  are  carried  out  for  ^(h)  as  well  as  for  A(h): 

A(h)  =  [  d£h  exp ( " ifl h  *  Qh>  A  (h)* 


a  -2 

(so  that  A(h)  =  (2rr) 


dQh  exp(i£h  •  Qh)  A(h)), 


cp(h) 


=  j"  dk^  exp  (ikj 


Qh)  *(h); 


A  r-' 

4(h)  -  ) 


Jn  +  V  <Vh’  Jv  <Vh>  »>-T>(-inh*h)r(h)  • 

h  h  h 


(7.13) 


where  ik  =  arctan  (q^ /q*),  and^'(h)  depends  on  v,  instead  of  n,.  Then 
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A(h)  =  £(h)  +  6  (G  .) 


,  oo  2tt 

2-  -  dq, 


^7^7T-Jd*u\(,!uahleJtp(-inh*u)Id^h 


9  j/_  mi.  v*  j.  nh  o 


[mwk  'lau^h1  +  ^  al-]  8(ah+ih  •  au>  »(Mi  (au) 


where 


(7.  14) 


00  00 


A^JsJ  dau  J  dpu  2TTmVauJ0(quau)  X'(u)  lv  .  0  , 

O  -00  u 


and  we  have  used  the  relation 


I  J„+v  J„  “  **'<*>■ 


(7.15) 


Furthermore,  eq.  (7. 13)  yields 

oo  oo 

I  (ah>  =2  I  da 

r 

"h 


h  J  dph  2nm2“2ai 

O  -00 


Jn  ^qhah^  exP(inhYh)  A  (h)- 
h 


(7.16) 


Therefore  eq.  (7.14)  may  be  reduced  by  means  of  this  barring  operation  to 
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CsJ 


2 

-  .  -  .  .  ,  4ne  N  \ 

”  a  %  inmuuA  L 


T  d\  I  dph  2M%Jn 

Jo  h  -oo  n  h 


exp  (in^  t^)  6_(Gx- [(kp^/m)  + 


rq  dq  2n 

- _H — H — -  r  J  (q  a 

q2  +  k2  +  h  o  U  "h  U 


exp  (-11^ *u)  [mu)k  ^  +  K%) 


*r5Tl{|V  h’VSu-V1*3/ 

h  h 


This  equation  can  be  written  in  a  notation  that  may  be  more  familiar 

A  (2)=  H(q)+J  J  dS'  K^.a')  A  (a')> 


-00  -00 


where 


K(a.a'>  =  7TXT  I  I  dpl 

iA(q  +k  +h  )  o  -® 


y  6_(G  -  [kp|/m  +  nu)]) 


n  =  -oo 


nl-^l  J  (J-JJ*)  exp  (in[*  -  *'] 
n  mu)  n  muj 
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SP  I! 


nmut' 


o 

?'P. 


+  mi’  (3Va)| 


■  S(£,k  -  3'  -  3),  (7-19) 

and  £  is  the  kinetic  momentum. 

Equation  (7.17)  or  (7.18)  is  of  the  Fiedholm  type,  but  the  kernel  K  is  rather 
complicated.  We  give  it  some  study  in  the  next  section,  but  no  detailed  examination 
has  yet  been  carried  out. 
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SECTION  8 

COMPARISON  OF  RESULTS  WITH  EXISTING 
KINETIC  EQUATIONS 

Our  results  reduce  to  the  homogeneous  kinetic  equation(1. 11. 17),  or  (1,11.1), 

(II.  4.  9),  when  cp(j)  is  independent  of  Qj.  The  easiest  way  to  show  this  is  to  return 
to  eq.  (7.1),  and  perform  the  integrations  over  the  j  Qj  |  .  We  obtain  factors 
6 ^  ),  6(^  -  £  ),  etc.  ,  so  that  the  integrations  over  all  but  one  of  the 

|  jf,  |  are  trivial;  we  may  set  =  ‘  •  •  =  ^ anc*  ^  =  ^u  =  "  ’  '  ^en  t*ie 

J 

R^sibois  summation  procedure  yields  geometric  series,  and  the  final  :  esult  is 
precisely  eq.  (II.  4.  91 

In  order  to  see  that  this  is  plausible  without  going  through  the  details,  one 
should  note  that  eq.  (7.17)  reduces  to 

Mah)=  s  (ah)/n  <*&•.=%> 

when  cp(h)  =  6(k,  )cp(h),  if  the  normalization  of  cp  is  changed  by  a  factor  £  .  Roughly 
speaking,  the  term  q(a)/r](a)  of  eq.  (II.  4.  9)  comes  from  the  term  A  (a,  h)  of  eq.  (7„  8), 
and  the  integrand  bJGaj[q^)(u)/r\{u)]/r}f  {u)  in  the  other  term  comes  from  the  factor 
r(u,  h;x)  in  eq.  (7.  8). 

Let  us  now  return  to  the  diffusion  problem.  So  far  as  we  are  aware,  no 
equivalent  results  have  yet  been  published.  The  only  comparable  results  are  those  of 
Eleonskii,  Zyryanov,  and  Silin^1^,  especially  their  equation  (17).  They  made  a 
quantum  mechanical  calculation  of  the  collision  integral,  using  the  firs*  Born  approxima¬ 
tion  and  the  Landau  representation  (Landau  levels),  and  later  took  the  classical  limit. 

The  collective  effects  were  taken  into  account  through  the  ad  hoc  introduction  of  a 
dielectric  tensor.  They  simplified  the  results  by  assuming  that  cp(j)  is  independent  of 

QX  ,  depending  only  on  Q^(or  y®  ,  in  our  previous  notation),  a  ,  p  ,  and  t.  If  we  set 
j  J  J  J  J 

the  collective  factor  equal  to  1  (so  that  the  results  are  analogous  to  the  Landau  equation), 

their  collision  integral  for  the  single  species  plasma  may  be  written  as  follows: 

oo 

WK  J 


2  2 

dp,  2tTm  uj  a, 
h  n 
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to  be  compared  with  our  result  (Cor  x -•  0); 
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cxp{i[n^  -nh][p. 


8h])exp(i[^  -  J J ]  Q, 


y 

;ah 


Expressions  (8.1)  and  (8.  2)  are  very  similar,  but  seem  to  be  slightly 

different.  We  have  not  yet  traced  the  origin  of  the  differences. 

2  2  ^ 

In  (8.1)  or  (8.2),  the  coefficient  of  3  /3Q^  ii  proportional  to  uu  ,  or 

-2  a 
B  .  Hence  the  results  may  be  in  qualitative  agreement  with  Fick's  law  with  a 

collisional  diffusion  coefficient.  (See  References  (1),  (2),  and  (12)  for  more  complete 

discussions.  )  Our  hope  is  that  by  properly  taking  into  account  the  collective  effects, 

as  we  have  done,  wo  may  eventually  be  able  to  understand  better  the  phenomenon  of 

“2 

anomalous  diffusion,  where  the  diffusion  coefficient  decreases  less  iapidlv  than  B 
for  increasing  magnetic  field  strength. 
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SECTION  9 
CONCLUSION 

We  have  provided  here  a  basis  for  kinetic  theoretical  studies  of  the  diffusion 
of  a  plasma  across  a  magnetic  field.  The  main  result  contained  in  this  report  is  a 
demonstration  that  the  problem  of  deriving  a  kinetic  equation  for  the  one  particle 
distribution  function  can  be  reduced  to  the  problem  of  solving  an  equation  of  the  Fredholm 
type  (7.18).  The  demonstration,  which  is  believed  to  be  original,  contains  no  explicit 
restriction  on  the  amplitude  or  length  scale  of  the  inhomogeneity  across  the  magnetic 
field. 

Much  further  work,  however,  remains  to  be  done  on  the  refinement  of  the 
physical  model,  and  on  the  qualitative  and  quantitative  examination  of  the  properties 
of  the  resulting  equations. 
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^Errata: 

Pane  1977,  line  26;  the  last  two  lines  of  the  expression  for  6L,  should  read 
,  n  hu 

as  follows: 

^[rru  H  Z3/6p^  -  (  £X  X  ?> / 3  x^  )/' 

+  |A(cos  [5^  -  6]  4/da^ 

-[l/a^.lsin[9^-'-t  ]d/39^J-[h  4 — >  u] 

Page  2000,  line  6  and  line  9:  "(8.15)"  should  be  "(8.14)"  . 
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f  Errata: 


Page  2008,  line  12:  "q  =  1,2,  2,  2,1"  should  be  "q  =  1,2,  2, 1,1".  Se  ■  also 
the  footnote  on  page  130  of  paper  II. 

Page  123,  line  5:  "m2"  -hould  be  "rr/". 

line  10:  "(2n)'1"  should  be  "(2tt)‘2". 

s  - 1 

Page  125,  Fig.  2:  "11"  should  be  "  f]  ". 

CT  n  =1 

Page  128,  line  6:  "m  "  should  be  "rn  >n. 

2  2 

Page  128,  last  paragraph:  "first  order  e  "  should  read  "first  order  in  e  ". 

Page  129,  line  21:  "  ^  ^(Gn^cpeC^(a)cpe<^(h)"  should  be  "  ^6(Ga^)Ga^cpe^(a)cpe<^(h) 

( v  V-1 

Page  129,  line  26:  "  — — -  "  should  be  " 
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ERRATA  (as  of  12  January  1966) : 


Page  2,  fourth  line  after  eq.  (2.1):  replace  "effect"  by  "affect". 

Fege  2,  third  paragraph:  replace  "assumption  (1.3.26)"  by  "assumptions  (1.3.26) 

and  (1,8.12)". 

Page  2,  end  of  third  paragraph:  replace  "with  W  set  equal  to  zero"  by  "with 
W  chosen  so  as  to  cancel  the  contribution  of  V(kz-0)  to  the  inter¬ 
action  matrix  elements." 

Page  2,  second  last  line  before  footnote:  replace  "may  be  a  less  trivial 
problem"  by  "is  a  less  trivial  problem". 

Page  7,  add  at  bottom  of  page: 

"It  should  be  observed  that  the  matrix  elements  that  are  diagonal 
with  respect  to  the  k's  are  not  necessarily  zero.  However,  in 
accordance  with  the  discussion  in  Section  2,  we  ignore  these  contributions 
by  replacing  (£^  +  fk^-k^]^  +  k^)  ^  by  0  when  k^  ■  ." 
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Page  8, 
Page  9, 
Page  11, 

Page  12, 
Page  14, 
Page  17, 


Page  21, 
Page  22, 
Page  22, 
Page  25, 
Page  25, 
Page  26, 


third  line  of  eq.  (4.1):  "{zj"  should  read  "{zj}". 

left  side  of  eq.  (4.5)  should  read  "f  (0)". 

r,s 

third  line  from  bottom:  replace  £  by  I  ,  and  include  a  marginal 

k  k(?0) 

note  referring  to  the  paragraph  added  at  the  bottom  of  page  7. 

line  3:  "k  ■  k"  should  read  "k  »•  -k". 
u  u 

second  line  of  eq.  (5.6):  replace  £  by  £ 

k  k(*0) 

add  a  dagger  after  “binary  correlation  function",  and  write  the  following 
foofnc-..  at  the  bottom  of  the  page: 

"tThe  appearance  of  at  this  stage,  instead  of  j  i  implies  that 

the  self-consistent  electric  field  is  being  ignored  (of.  Sections  2  and 
3)." 

line  2:  "G  "  should  read  "G 
s  o 

eq.  (7.8):  replace  "dG^"  by  "(dGx/27r)". 

eq.  (7. 1C):  replace  "$(h)  Ma,h;x)"  by  "$(h)  A(a,u;x)". 

eq.  (7.18):  "¥($)"  should  read  "5(g)". 

add  "  )  "  at  end  of  last  line. 

last  line  of  eq.  (7.19):  replace  "k"  by  "k^'  . 
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Pa^e  28, 

Page  29, 
Page  29, 


eq.  (8.1):  replace 
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\<Vh> 


2  by  |jdty(k2  +  t2)"1 


exF(l[na  -  njb)  expUt*^)  J  ».„)  J  «VU 

a  n 

line  3:  "k3/3p  "  should  be  "k3/3p 


lines  3  and  6:  replace  "Expressions  (8.1)  and  (8.2)  ...  differences." 
by  "The  results  (8.1)  and  (8.2)  are  identical." 


Page  31,  add  at  bottom  of  page: 

"Page  1996,  eq.  (7.21):  replace  Z(e^)  by  Z(0),  and  add  '+  0(eSv. 
Eq.  (7.22) :  add  '+  0(eA)'. 
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Lines  15  to  17:  replace  'and,  as  it  happens,  is  Independent  of  e  7 
by  'to  first  order'." 

Page  32,  last  correction:  replace  "d"  by  "D"  in  two  places. 
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